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Abstract 

Context. We seek clarification of the nature of X-ray sources detected in M 31. Here we focus on CXOM31 J004253. 1+41 1422, the 
brightness of which suggests that it belongs to the class of ultraluminous X-ray sources. 

Aims. We determine the X-ray properties of sources detected in the XMM-Newton / Chandra monitoring program. We investigate 
spectral properties and search for periodic or quasi-periodic oscillations. A multi-component model is applied to the spectra obtained 
from XMM-Newton data to evaluate the relative contributions from thermal and non-thermal emission. The time dependence of this 
ratio is evaluated over a period of forty days. 

Methods. We simultaneously fit data from XMM-Newton EPIC-pn, MOS1 and MOS2 detectors with (non-thermal) powerlaw and 
(thermal) multicolored blackbody. 

Results. The X-ray spectrum is best fit by the combination of a thermal component with kT ~ 1 keV and a powerlaw component with 
photon index approximately 2.6. From combined analysis of Chandra, Swift and XMM-Newton data, the unabsorbed total luminosity 
of this source decreases from ~ 3.8 x 10 3 ' erg s in the first observation to ~ 0.5 x 10 39 ergs s over a period of three months. The 
decay closely follows an exponential decline with a time constant of 32 days. The source spectrum evolves significantly, exhibiting a 
faster decline of the thermal component. We do not find evidence of any significant temporal features in the power density spectrum. 
The presence of a thermal component at kT ~ 1 keV in conjunction with a non-thermal high energy tail, is also consistent with 
spectral properties of other ULXs in the "high state". 

Conclusions. Our analysis indicates that the underlying source of this first ULX in M 31 is a black hole of mass, M > 13 Mq, 
accreting near the Eddington limit, that underwent a transient outburst followed by an exponential decay reminiscent of transients 
associated with galactic X-ray novae. 
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1. Introduction 

Ultraluminous X-ray sources (ULXs) are very bright point 
sources with an X-ray luminosity of L t > 10 39 erg s _1 , which 
exceeds the Eddington luminosity for compact objects with mass 
of approximately 10 Mq. This source class was first identified by 
the Einstein Observatory in the eighties (Fabbian o|1989 1. These 
sources are not associated with the center of galaxies and thus 
do not belong to the class of super-massive black holes. On 
the other hand, they are too bright to be associated with stel- 
lar mass black holes with M < 10 Mq, accreting below the 
Eddington rate. An exciting possibility is that the underlying 
sources of ULXs are intermediate mass black holes (IMBHs), 
in the mass range 10 2 - 10 4 Mq, accreting at sub-Eddington 
rates (Mi ller & Colbert(2 004). For example, observations of an 
ULX in NGC5408 suggested a mass of ~ 100 Mq, indicating a 
possible existence of IMBHs (Fabian & Ward 1993| l. However, a 
recent study of this source's variability indicates that the mass is 
likely below 100 M Q ( |Middleton et al.|201 l( l. Although the iden- 
tification of ULXs as IMBHs has not yet been firmly established, 
one should contemplate the theoretical progenitors, that could 
result in producing such objects. Several theories have been sug- 



* Based on observations obtained with XMM-Newton, an ESA 
Science Mission with instruments and contributions directly funded by 
ESA Member States and NASA. 



gested for the formation of IMBH. |Madau & Rees| p001 ) stud- 
ied the collapse of Population III stars, while others considered 
the collapse of massive stars (in young super-massive star clus- 
ters) e.g. [Portegies Zwart et al. ( |2004| l. Another possibility is 
of course that the masses are much less than 100 Mq, but that 
the source accretes at Super-Eddington rate, as suggested by 
Begelman (2002). Yet another factor to be considered is devia- 
tions from the common assumption of isotropic emission as sug- 
gested by King et al. ( 200 1 )l . A long these lines for anisotropic 



emission, Reynolds et al. ( 1997 1 proposed a link between ULXs 
and extreme beaming, as observed in galactic micro-quasars. 
The nature of the ULX class remains unclear, which is in part 
due to the fact that the sample size is still rather small. 

To better understand the ULX phenomenon, an increase in 
sample size and careful consideration of spectral and sample 
properties is desirable. The ULX sources are mainly found in 
star forming galaxies. Until recently, no ULX had been identified 
in M 31. Our ongoing monitoring program^] for resolving super 
soft source states of optical novae in the central area of M 31 
(P.I. Wolfgang Pietsch) with XMM-Newton and Chandra HRC- 
I yielded a transient source, which, as discussed in this paper, 
represents the first recognized ULX in our companion galaxy. 



www.mpe.mpg.de/m31novae/xray/index.php 



1 



A. Kaur et al.: CXOM31 J004253.1+411422:The first ULX in M 31 



We present the time development of the outburst and the 
spectral analysis of a ULX in M 31 using XMM-Newton, 
Chandra and Swift data over a period of approximately three 
months (Sect. 2), followed by the summary of the results in Sect. 
3. Conclusions are presented in Sect. 4. 




Figure 1. XMM-Newton EPIC image of the central part of M31 
produced by combining pn, MOS1 and MOS2 data from all 5 
observations. Red, green and blue show the (0.2 - 1.0) keV, (1.0 
- 2.0) keV) and (2.0 - 12.0) keV energy bands, respectively. 



2. Observations and data analysis 

CXOM31 J004253. 1+41 1422 was discovered with 
Chandra HRC-I in a 19 ks observation on 2009 Dec 17.89 UT 
(Henze et al. 2009). The source position was determined to be 
R.A (J2000) = 00:42:53.15, Dec. (J2000) = +41:14:22. 9 using 
the catalog of X-ray sources in M 3 1 assembled by |Kaaret 
2002) for relative astrometry. Figure 1 shows the position of 



the source with respect to other X-ray sources near the center 
of M31 in an XMM-Newton color image. A transient optical 
counterpart with m(F435W) = 23.8 mag was identified in 
HST observations ( |Garcia et al.||2010| . The X-ray source was 
first detected with Chandra HRC-I. Level 2 event files were 
analyzed using (Chan dra Interactive Analysis of Observations; 
Fruscione et al.|2006| P1 to obtain count rates as shown in Table 
1. An adapted version of the XMMSAS tool emldetect was 
used to estimate background-corrected and exposure-corrected 
fluxes and count rates (see |Henze et al.||2010j i. Continuous 
monitoring of this source was then carried out by Swift-XRT 
from December 22 - 27, 2009. The data were analysed using 
the HEAsoft XIMAGE package (version 4.5.1) with the sosta 
command (source statistics) for estimations of count rates. 
We took into account the XRT PSF of the sources that we 
computed with the command psf, as well as exposure maps 
that were created with the XRT software task xrtexpomap 
within XIMAGE. The count rates obtained from Swift-XRT 
and Chandra HRC-I before and after the XMM-Newton 
observations were converted to unabsorbed fluxes using energy 
conversion factors (ecfs), which were computed using fakeit 
in XSPEC assuming the best fitting spectral model for the 
first and last XMM-Newton observations, respectively, and 

using publicly available instrument response files. We obtain 

ec f . . — A, 1 v 1f\l0 „ tt , a ,.„-l „„-,2 



lHRC-I-1 



6.3 x 10 cts erg cm and ecf hrc-i-2 = 7.7 x 



10 10 cts erg 1 cm 2 , respectively. For the Swift-XRT data, the 
ecf has been computed to be 1.57 x 10 10 cts erg -1 cm 2 . 



The spectroscopic data were then obtained by 
XMM-Newton from December 28, 2009 until February 
02, 2010 in five distinct observations using the European Photon 
Imaging Ca mera (EPIC) EPIC -pn ( |Striider et al.|2001| l, MOS1 
and MOS2 ( |Turner et al.|2001| ) CCD detectors are mounted on 
the three X-ray telescopes on XMM-Newton. The exposure 
times as well as the raw count rates obtained with these detectors 
are presented in Table 1. The XMMSAS version 10.0.0 was used 
to filter the standard pipeline event files, to generate images, 
light curves, spectra and the appropriate detector response 
functions. The event files for all three cameras, pn, MOS1 and 
MOS2 were analyzed with the following parameter settings: 
For spectrum generation, we set FLAG = to reject bad pixels, 
and PATTERN < 4forpnand< 12 for MOS1 and MOS2 event 
files were allowed to reduce noise in the data. Light curves were 
derived to investigate the variability in the source. Barycenter 
correction was performed using the barycen task in XMMSAS. 
The light curves were then corrected from various effects such as 
bad pixels, GTIs, vignetting using epiclccorr. As an example 
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Figure 2. Light curve from XMM-Newton EPIC-pn for obser- 
vation 0600220201 with time binning of 100 s in 0.2- 10 keV 
energy band. 

Fig. 2 displays the observed count rate for XMM-Newton 
EPIC-pn for the first observation, binned in 100 s intervals. The 
power density spectra were generated using the task powspec 
in FT00LS using the Fast Fourier transform algorithm (FFT). 
The normalization was chosen such that the white noise level 
expected from the data errors corresponds to a power of 2. For 
spectrum extraction, the XMMSAS tool epatplot was used to 
check for pile-up for appropriate source extraction region. The 
inner PSF part of these regions were excluded to avoid pile-up. 
Response files were generated by using SAS tools rmfgen 
and arfgen. The spectral binning was constructed to attain 20 
counts per bin to assure uniform statistics across the energy 
range. Bad pixels were ignored and energy channels from 0.2 
- 10.0 keV were considered for spectral fitting using XSPEC 
version 12.6.0. 

3. Results 

Our initial spectral fitting for XMM-Newton data used a bro- 
ken powerlaw model (BKNPOWER in XSPEC), but the ^ 2 /d.o.f. 
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Table 1. Observations log 



Telescope/instrument 


Obs ID 


Date 
(UT) 






Exptime 
(ks) 


Rate" 
(ct s" 1 ) 




L* 

ninth 

(1CF 


erg s 1 ) 


Chandra HRC-I 


10885 


2009- 


•12- 


08.94 


18.27 


< 1.5e-03 


0.002'' 


Chandra HRC-I 


10886 


2009- 


12- 


•17.90 


18.34 


3.300 ± 


0.040 


3.77 


±0.04 


Swift-XRT 


00031518013 


2009- 


12- 


•22.04 


3.6 


0.660 ± 


0.200 


3.04 


±0.20 


Swift-XKT 


00035336016 


2009- 


•12- 


23.05 


4.2 


0.720 ± 


0.100 


3.31 


±0.10 


Swift-XRT 


00035336017 


2009- 


•12- 


24.04 


4.8 


0.620 ± 


0.100 


2.86 


±0.10 


Swift-XRT 


00035336018 


2009- 


-12- 


25.05 


5.2 


0.600 ± 


0.100 


2.76 


±0.10 


Swift-XRT 


00035336019 


2009- 


12- 


26.26 


5.0 


0.560 ± 


0.100 


2.58 


±0.10 


Swift-XRT 


00035336020 


2009- 


•12- 


•27.07 


5.0 


0.580 ± 


0.100 


2.67 


±0.10 


XMM-Newton 


0600660201 


2009- 


12- 


•28.53 


16.86 


6.447 ± 


0.022* 


2.16 


± 0.07^ 














1.810 ± 


0.010 c 


















1.851 ± 


0.0 10 rf 






XMM-Newton 


0600660301 


2010- 


01- 


07.32 


15.43 


1.784 ± 


0.012 


1.49 


± Q.Q5 f 














1.105 ± 


0.008 


















1.143 ± 


0.008 






XMM-Newton 


0600660401 


2010- 


01- 


15.53 


15.33 


3.832 ± 


0.017 


1.16 


± 0.04 f 














1.372 ± 


0.009 


















1.402 ± 


0.009 






XMM-Newton 


0600660501 


2010- 


01- 


25.11 


17.83 


3.042 ± 


0.015 


0.71 


± 0.03' 














0.846 ± 


0.007 


















0.873 ± 


0.007 






XMM-Newton 


0600660601 


2010- 


•02- 


02.11 


15.43 


1.072 ± 


0.009 


0.65 


± 0.03' 














0.694 ± 


0.006 


















0.687 ± 


0.006 






Chandra HRC-I 


10808 


2010- 


•02- 


•15.86 


17.12 


0.547 ± 


0.080 


0.51 


±0.08 


Chandra HRC-I 


11809 


2010- 


•02- 


26.27 


18.42 


0.490 ± 


0.180 


0.46 


±0.18 



Notes. <n) Count rate as observed from the source without rejecting pile-up pixels. {b) EPIC-pn, (c) MOS1, {d) MOS2 data. 3cr upper limit 
Derived from the simultaneous fitting of EPIC-pn, MOS 1 and MOS2 in XSPEC. te) Total unabsorbed luminosity from the source within 
(0.2-10 keV) energy band. 

Table 2. The spectral parameters obtained from the model fitting from XMM-Newton 



OB SID 


N" 
(10 20 


31 

cm" 2 ) 


kT b 
(keV) 




P 




* 2 /d.o.f. 


Rf„ V«w(i) 
(km) 


(10 39 erg s _1 ) 


(10 39 erj 


BB 

is' 1 ) 


(10 39 erg s-') 


0600660201 


5.1 ± 


1.1 


1.070 ± 


0.010 


2.59 ± 


0.16 


1318.87/1165 


52.90 ± 0.01 


0.33 ± 0.07 


0.94 ± 


0.07 


2.16 ± 0.07 


0600660301 


7.1 ± 


1.2 


0.993 ± 


0.016 


2.58 ± 


0.13 


963.62/991 


58.89 ± 0.02 


0.32 ± 0.05 


0.71 ± 


0.05 


1.49 ± 0.05 


0600660401 


9.1 ± 


1.2 


0.913 ± 


0.013 


2.81 ± 


0.14 


1046.74/953 


57.32 ± 0.02 


0.33 ± 0.04 


0.36 ± 


0.04 


1.16 ±0.04 


0600660501 


5.3 ± 


1.0 


0.815 ± 


0.014 


2.51 ± 


0.12 


962.61/938 


61.57 ±0.02 


0.23 ± 0.03 


0.27 ± 


0.03 


0.71 ± 0.03 


0600660601 


6.5 ± 


1.6 


0.769 ± 


0.013 


2.70 ± 


0.13 


830.66/821 


64.18 ±0.03 


0.16 ±0.03 


0.25 ± 


0.03 


0.65 ± 0.03 



Notes. ((!) External absorption column density. (w Temperature of inner-disc from DISKBB. (r) Photon index from powerlaw model. << " R,„ is inner 
radius of the multicolored blackbody disk, and i is the inclination angle. <f) Absorbed luminosity contribution from the powerlaw model (hard 
component) in (0.2-10 keV) energy band. Absorbed luminosity from the multicolored disk component within (0.2-10 keV) energy band. 
(s) Total unabsorbed luminosity from the source within (0.2-10 keV) energy band. 



= 2953.96/1174 clearly indicated that this is not a satisfactory 
model. The luminosity derived from this simple fit indicated the 
fact that this source is a member of the ULX class. For these ob- 
jects, an alternative, often better fitting model is a combination of 



a non-thermal (power law) and thermal component (Gladstone 
|et al.|2009] l. To explore this possibility we applied a model that 
combines a non-thermal powerlaw component (PO) with a mul- 
ticolored blackbody component (DISKBB). The thermal compo- 
nent is associated with an accretion disk around a central black 



hole as developed by |Mitsuda et al. ( 1984 1 and Makishima et al. 
(1986). This model gave a significantly improved fit, indicated 
by x z I d.o.f. values in Table 2. 



In addition to the intrinsic two component model, the source 
flux must be corrected for extinction by foreground gas in the 
Milky Way and the local gas along the line of sight through 
M 31. Extinction was modeled using the Tuebingen-Boulder 
ISM absorption model (TBABS) ( | Wilms et al.|2000| l. The hydro- 
gen column d ensity ascribed t o the Milky Way was fixed to 5.32 



x 10 cm (Dickey & Lockman 1990 1, whereas the column 



density for the host galaxy was included as a free parameter. 
However, the derived value of the column density in M 3 1 must 
thus be considered with caution because we do not have informa- 
tion about the actual metallicity along the line of sight through 
M 31 and simply assume the Milky Way template expressed in 
terms of a standard cross-section per hydrogen atom. 
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To illustrate this fitting procedure, Fig. 3 shows the 
XMM-Newton EPIC data for first of our five observational 
epochs. The decomposition of the X-ray spectrum into the 
model components in physical units is shown in Fig. 4. 




1 2 
Energy (keV) 

Figure 3. Combined powerlaw and multicolored blackbody fit 
to the joint set of data from EPIC-pn (black), MOS1 (green), 
MOS2 (red) with X 2 I d.o.f. = 1318.87/1165 for the observa- 
tion 0600660201 with XMM-Newton. The model parameters 
obtained from this spectral fit are presented in Table 2. 



and Swift-XRT. The total unabsorbed luminosity declined from 
about 3.7 x 10 39 to 0.5 x 10 39 erg s~' during the period of ob- 
servation. We note that this source was initially more luminous 
than the total X-ray lumi nosity of M 31 in t he (0.1-2.4) keV 
band sampled by ROSAT (Suppei ^et al.|1997) . Futhermore, this 
source was atleast one order of magnitude brighter than any 
of the 45 X-ray transients detected in M 31 by Chandra and 
XMM-New ton from October, 1999 to August, 2002 ( |Williams 
|et al.|2006| ). 
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Figure 4. The theoretical model underlying the data fit (solid 
line) in Fig. 3 is shown in decomposition between the thermal 
(dash-dotted line)and the non-thermal (dotted line) component. 



To convert fluxes to luminosities, we assume a source dis- 
tance of 780 kpc ( |Holland|[T998l |Stanek & Garnavich|[T99"8] l 
and isotropic emission. The resulting luminosities during the 
XMM-Newton observations are about 10 39 erg s _1 , which 



places this source in the mid range of known ULXs (Gladstone 
|et al.|2009) . Table 2 summarizes the derived luminosities (from 
XMM-Newton only) as well as the effective inner disk radius. 
As a function of time, both the non-thermal (Lpo) and ther- 
mal (L BB ) component decrease, but the latter decreases more 
rapidly (see Fig. 5). Table 1 lists the calculated luminosities de- 
rived from count rates associated with observations of Chandra 



Figure 5. Temporal variation of the total unabsorbed luminosity 
(Lunabs) as per observations of XMM-Newton. The decomposi- 
tion into the powerlaw (L po ) and blackbody component (L^) re- 
veals a distinct evolutionary behavior of these two components. 
The time zero corresponds to Chandra HRC-I observation ID 
10886 (see Table 1). 

From the X-ray spectrum fits of the first XMM-Newton ob- 
servation (when the source was nearly brightest and we had the 
best quality data), we extract an estimate of the mass of the black 
hole following the formalism given by Makishima et al. (2000). 
Under the assumption of a geometrically flat and optically thick 
accretion disk, these authors derived the following relation: 



M = 



R;, 



8.86 a 



M Q , 



where R,„ is the inner radius of the accretion disk, measured 
in km and a is a dimensionless parameter that relates the inner 
radius to the Schwarzschild radius via R,„=3a' Rs- 

For the bolometric luminosity emitted by the accretion disk, 
Makishima et al.|(|20Q0l> find (their eqt. 9): 



( 



M 



\l0M o ) \keVl 
ergs~ 



where T,„ is the temperature characterizing the inner accre- 
tion disk. The remaining parameters specific to their model are 
here chosen to take their standard values as indicated by the nor- 
malization. The spectral fit yields the products of R,„ and the 
square root of the cosine of the inclination angle of the disk, as 
well as temperature, T,„, and the bolometric luminosity (which 
we assume to be unabsorbed L/,/, component (= 1.6 x 10 39 erg 
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s ) integrated from 0.2-10.0 keV). From the two equations 
above, we can therefore determine the mass and inclination an- 
gle, although significant uncertainties in both quantities are as- 
sociated with the model parameters, a, £ and k fixed to the val- 
ues 1, 0.41 and 1.7, respectively (see Makishima et al. 2000). 
The a parameter has been defined above and £ is a correction 
factor, reflecting the fact that T,„ occurs at a radius somewhat 
larger than R,„. The specific value adopted here was motivated 
by Kubota ~et al.| ( |1998] l. The final paramter, k represents the ra- 
tio of color temperature to effective temperature and the value 
adopted here has been taken from Shimura & Takah a"ra| ( |1995| ). 
With these caveats in mind, we find a mass of 13 M0 with a 
statistical uncertainly of about 4 % . This implies an inner ra- 
dius of R,„ = 124 km. From the fitting parameter R,„ ~\Jcos(i) (see 
Table 2), one then infers an inclination angle of i ~ 80°. The as- 
sumption of a non-rotating Schwarzschild black hole (a = 1) 
can lead to a significantly underestimated black hole mass as the 
extreme case of a maximally rotating Kerr black hole {a - 1/6) 
allows for a more massive central object approaching the realm 
of IMBH. Note, that this increase with black hole spin relies on 
the assumption that the inner edge of the accretion disk can be 
identified with the location of the innermost stable circular orbit 
(ISCO). However, determining the inner edge of a black hole ac- 
cretion disk is a far more complex issue, e.g 
( |20T0l >. 

As described in Sect. 2, we searched for an underlying peri- 
odicity in each of the observing runs over a range of frequencies 
limited by the temporal extent of the observation and the bin- 
ning time of 100 s. The observed count rates (e.g. Fig. 2) display 
fluctuations within the range expected from counting statistics. 
None of the power density spectra reveal a significant feature. 



We compared various properties of this ULX with a se- 
lected sample from the XMM-Newton archival study of ULXs 
by Winter et al. (2006) in 32 nearby galaxies with distance < 



Abramowicz et al. 
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Figure 6. Temperature vs Luminosity for 19 selected ULXs from 



Wint er et aLl (f2C)06 ) along with the candidate ULX (red) in M 3 1 . 
The M 31 source is a member of the "high temperature" sub- 
group. 



4. Discussion and conclusions 

The intrinsic luminosity of the source is consistent with those 
found for other ULXs in nearby galaxies. The presence of a ther- 
mal component at kT ~ 1 keV in conjunction with a non-thermal 
high-energy tail, is also consistent with spectral properties of 
other ULXs. The temperature and photon index of M 31 source 
fall in the center of observed distributions for ULXs ( |Winter| 
et al. 2006), lending additional support to a ULX nature. 



8 Mpc and unabsorbed luminosities, L x > 10 erg s in the 0.3 
- 10 keV energy range. Motivated by the observed properties 
of black hole binaries in the Milky Way, these authors classi- 
fied their sample into two subsets, "low state" and "high state", 
based on their respective spectral properties. Low state sources 
were characterized by a single powerlaw X-ray spectrum, while 
high state sources required an additional multicolored blackbody 
with powerlaw. Since the spectra described in Sect. 3 clearly ex- 
hibit the best fit when a powerlaw is combined with multicolored 
blackbody, this source should be classified as "high state" ULX. 
Therefore we limit our comparison to 19 of high state ULXs 
in their sample. This sample is slightly reduced because we se- 
lected the brightest epoch out of multiple XMM-Newton ob- 
servations, and those were not always available in [Winter et al. 
( |2006> . 

Figure 6 displays the resulting set of "high state" ULXs in 
the L-T plane. As Winter et al. (2006) pointed out, the tempera- 
ture distribution is possibly bimodal. The ULX in M 3 1 belongs 
to the high-T group, which is commonly believed to represent 
the class of stellar mass black holes. The low-T group is char- 
acterized by a larger range of X-ray luminosities and has been 
tentatively attributed to the hypothetical class of IMBHs ( | Winter 
|et al.||2006| . Although the sample size is still very small, it ap- 
pears that a physical gap between these two groups exists. If that 
is the case, an explanation may be related to different formation 
scenarios for these groups. As Winter et al.| ([2006 ) suggested, 
one possibility is associated with the unique features of pair in- 
stability supernovae (PISNe) in the early universe (Population III 
origin, Heger et al. 2003). The mass of ~ 13 Mq derived in Sect. 
3 suggests that we are dealing with a normal stellar remnant, one 
that could have been formed later in the universe, provided a is 
not significantly less than unity due to rotation. 



> 

CD 



o 



CM 
O 




Figure 7. For the same sample shown in Fig. 6, a comparison of 
photon index vs temperature reveals that the M 31 ULX cools 
significantly, but remains within the range covered by the "high 
temperature" group. The labels 1 through 5 correspond to the 
chronological entry of XMM-Newton observations as listed in 
Table 1. 



Figure 7 shows the "high state" ULXs distribution in the 
T— r plane. The low-T group shows a large range in photon in- 
dex, resembling the large range in luminosity, however we note 
that the uncertainties in this parameter are very large. In con- 
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help with scheduling the TOO observations. M. Henze acknowledges support 
from the BMWI/DLR, FKZ 50 OR 1010. 
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Figure 8. Variation of the total unabsorbed luminosity (L una b s ) of 
ULX in M 31 as observed by Swift-XRT, Chandra HRC-I and 
XMM-Newton resembling an exponential decay with a time 
scale of 32 days. The time zero corresponds to Chandra HRC - I 
observation ID 10886 (see table 1). 



trast to the previous figure, here we show the temporal evolu- 
tion of the M 3 1 ULX, while for other ULXs, no evolution has 
been shown. Over a time span of forty days of observations with 
XMM-Newton, the spectral evolution clearly exhibits cooling, 
but almost no change in the shape of the power-law component. 
Referring back to Fig. 5, the drop in the temperature is also ac- 
companied by a significant change in the luminosity of this com- 
ponent. The power-law component, on the other hand, exhibits 
a declining luminosity as well, but not as rapid as the thermal 
one. The overall evolution of this transient is thus a combination 
of a strongly varying thermal component and a significantly less 
varying non-thermal part. 

Combining all data obtained with Chandra, Swift and 
XMM-Newton indicates that the X-ray luminosity of this 
source follows closely an exponential decline (see Fig. 8) with 
a time constant of 32 days. A continuation of this trend is evi- 
dent from data taken 150 days afte r the outburst when the source 
was detected at ~ 10 38 erg s" 1 ( Barnard et al.|[20lT l. This ob- 
served exponential decline of the luminosity is consistent with 
those determined for the class of galactic X-ray novae (Chen 
et al.|199~ I, which show FRED (Fast Rise Exponential Decay) 
like light curves, though the FR part in our case was missed. 
However, the light curves of this source class often exhibit more 
complex behaviour than a simple exponential decay. 

In summary, the Chandra discovery of 
CXOM31 J004253. 1+41 1422 has established it as the first 
member of the class of ULXs in M 3 1 . Follow up observations 
with Swift and XMM-Newton revealed an exponential decline, 
reminiscent of the late time evolution of Galactic X-ray novae. 
Spectral analysis of the XMM-Newton data suggests that this 
source was in a "high state" at the time of observation, and 
that the underlying source is likely a stellar mass black hole 
accreting near the Eddington limit. 
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